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Ph . Summary. We study aspects of the hydrodynamics of one-dimensional totally asym- 

^ ! metric iiT-exclusion, building on the hydrodynamic limit of Seppalainen (1999). We 

'^ I prove that the weak solution chosen by the particle system is the unique one with 

J3 ' maximal current past any fixed location. A uniqueness result is needed because we 

can prove neither differentiability nor strict concavity of the flux function, so we 

cannot use the Lax-Oleinik formula or jump conditions to define entropy solutions. 
K^ , Next we prove laws of large numbers for a second class particle in i^-exclusion. The 

(^ I macroscopic trajectories of second class particles are characteristics and shocks of 

^ ' the conservation law for the particle density. In particular, we extend to -Pf -exclusion 

Q^ ■ Ferrari's result that the second class particle follows a macroscopic shock in the Rie- 

O . mann solution. The technical novelty of the proofs is a variational representation 

for the position of a second class particle, in the context of the variational coupling 

method. 
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1. Introduction 

This paper continues the study of the hydrodynamics of totally asymmetric 
nearest- neighbor X-exclusion processes. These are interacting systems of indis- 
tinguishable particles on the one-dimensional integer lattice Z, subject to the re- 
striction that each site contain at most K particles. Independently at each site and 
at constant exponential rate one, one particle jumps to the next site on the right, 
provided the iiT-exclusion rule is not violated. 

A hydrodynamic limit for this process was proved in Seppalainen (1999). This is 
a law of large numbers according to which the empirical particle density converges 
to a weak solution of a nonlinear conservation law, under the appropriate scaling 
of space and time. 

The present paper begins with the question of uniqueness of the weak solution 
chosen by the particle system. The criterion we give is maximization of current 
across any fixed position. This issue arises because we do not know whether the flux 
function of the conservation law is strictly concave or differentiable. Consequently 
we cannot use the Lax-Oleinik formula to deflne an entropy solution, or one-sided 
bounds on jumps as an entropy criterion. 

We address this difficulty because it is closely tied to the main topic of the 
paper, the dynamics of a second class particle in the hydrodynamic scaling. As we 
shall see, our incomplete understanding of the macroscopic equation constrains the 
macroscopic description of the second class particle. 

A second class particle is an extra particle in the system that yields right of way 
to the regular, first class particles. The second class particle can jump only if there 
is no first class particle at its site who can jump, and a second class particle has to 
jump backward if first class particles fill its site. 

Second class particles in asymmetric exclusion processes have been studied for 
several reasons. In a certain technically precise sense, the second class particle 
marks the location of a shock at the microscopic particle level. These ideas were 
developed by Ferrari, Kipnis and Saada (1991), Ferrari (1992), and Ferrari and 
Fontes (1994a). At this same time, Ferrari (1992) discovered that in the hydrody- 
namic limit a second class particle follows the characteristics of the macroscopic 
equation. The connection between second class particles and characteristics was 
further elucidated by Rezakhanlou (1995) and Ferrari and Kipnis (1995). The 
second class particle also enters naturally in the study of the fiuctuations of the 
current, as in Ferrari and Fontes (1994b). 

Another strain of work, technically related to the investigation of microscopic 
shocks, concerns the invariant distributions of the exclusion process as seen from 
a second class particle, and the invariant distributions of the two-species process 
that has positive densities of both first and second class particles. This program 
culminated in the complete description of these measures by Ferrari, Fontes and 
Kohayakawa (1994). These and other results, such as the central limit behavior 
of the second class particle, are discussed in Part III of Liggett (1999) and in the 
notes of Chapter 9 in Kipnis and Landim (1999). 

The above results are for systems whose invariant distributions are explicitly 
known, namely the exclusion and zero-range processes. And the proofs typically 
use precise and deep knowledge about those invariant measures. 

In the present paper we prove laws of large numbers for a single second class par- 



information is currently available about the invariant distributions of these asym- 
metric X-exclusion processes. 

Of the papers mentioned above, Ferrari (1992) and Rezakhanlou (1995) are the 
closest predecessors of our paper. Both derived laws of large numbers for the 
second class particle. The limiting trajectory is a characteristic or a shock of the 
conservation law of the particle density. Ferrari proved strong laws in an asymmetric 
simple exclusion process, mainly in these settings: the initial particle distribution is 
either an i.i.d. product measure, or a product measure that represents a shock profile 
that is constant on either side of the origin (this is known as the Riemann problem). 
Ferrari's proofs rely on couplings of processes of several species of particles, and 
knowledge of invariant distributions. 

Rezakhanlou (1995) proved L^-laws of large numbers for the second class particle 
in asymmetric simple exclusion and in a totally asymmetric zero-range process. He 
assumed that initial distributions are product measures, and the initial macroscopic 
profiles are bounded integrable functions on R, with some additional technical 
assumptions. His proof derives the limit of a single second class particle from 
the limit of a macroscopically visible interval of second class particles, and the 
latter he gets from the basic hydrodynamic limits for exclusion and zero-range 
processes, from Rezakhanlou (1991). The invariant distributions of the processes 
enter his proofs through the macroscopic equations. These are known explicitly 
precisely because the invariant distributions are known. Rezakhanlou's approach 
would prove some result for i^-exclusion, since one can use the hydrodynamic limit 
of Seppalainen (1999) instead of Rezakhanlou (1991). 

Our approach is different from the earlier ones. We use a representation of 
the second class particle in the variational coupling method. To preview this, let 
us switch from the exclusion process to the associated marching soldiers process 
(in other words, consider the interface process instead of the increments process). 
In this approach a process z{i, t) from a general initial interface is coupled with 
auxiliary processes w^{ij t) that evolve from translated wedge-shaped initial profiles. 
The coupling is arranged so that z{i,t) is the envelope of the auxiliary processes: 

(1.1) z{i,t) = supw''{i-k,t). 

kez 

Let X{t) denote the position of a second class particle added to the exclusion process 
r]{i, t) = z{i, t) — z{i — 1, t). The new tool is this variational expression: 

(1.2) X{t) = inf{z : z{i, t) = w^{i- /c, t) for some k > X{0)}. 

Variational formula (1.1) is a microscopic analogue of the Hopf-Lax formula of 
viscosity solutions of Hamilton-Jacobi equations. In this same vein, (1.2) corre- 
sponds to the definition of forward characteristics through the Hopf-Lax formula. 
This justifies regarding the second class particle as a microscopic analogue of the 
characteristic, as expressed in the title of the paper. 

Variational formulas (1.1)-(1.2) permit us to work directly on the paths of the 
process, and we do not need a particular form for the initial particle distributions. 
We need to assume that the initial particle profile and the initial position of the 
second class particle both converge under the appropriate space scaling. The result 
is a law of large numbers for the second class particle. The theorems give either 



initially. The limit is a characteristic or a shock of the macrosopic partial differential 
equation. 

Here we run into the difficulty mentioned above. For 2 < i^ < oo (cases other 
than exclusion or zero-range) we do not have sufficient control of the regularity of 
the flux / of the conservation law. For this reason we cannot assert that macro- 
scopically the second class particle follows the solution of i; = /'(p(x,t)), where 
p{x,t) is the macroscopic particle density. Instead, we use the Hopf-Lax formula 
to define the macroscopic trajectory of the second class particle. 

If the macroscopic characteristics are unique (there are no rarefaction fans) we 
do not need the assumption of an initial law of large numbers for the second class 
particle. In this case the second class particle follows, in the hydrodynamic scale, 
the macroscopically defined characteristic that emanates from its random initial 
position. 

We improve earlier results in two ways: by broadening the class of processes be- 
yond exclusion and zero-range, and by permitting more general initial distributions. 
In one important sense the results of Ferrari and Rezakhanlou remain more gen- 
eral than ours: their results cover an asymmetric simple exclusion process where 
particles can jump both right and left, while our results are restricted to totally 
asymmetric processes where particles jump only to the right. This is because a way 
to apply the variational coupling approach without total asymmetry has yet to be 
discovered. 

Overview of the paper. In Section 2 we discuss the results. Section 2.1 reviews 
the hydrodynamics of totally asymmetric X-exclusion and zero-range processes 
(Theorem 1). We state a uniqueness criterion for the weak solution chosen by the 
particle system, valid for a concave flux that is not necessarily differentiable or 
strictly concave (Theorem 2). Section 2.2 states the laws of large numbers for the 
second class particle (Theorems 3 and 4). Section 2.3 comments on the various 
assumptions and conclusions. 

In Section 2.4 we discuss a generalization of the i^-exclusion and zero-range 
processes to which our results apply. In this model particles may jump in batches 
of varying sizes from one site to the next. Such a generalization is natural from a 
queueing perspective, where the particles are customers moving through a sequence 
of service stations. 

Some standard tools of conservation laws are not at our disposal for lack of 
information about the regularity of the flux / of the conservation law. This problem 
is clarified in Section 3. Without strict concavity of / the Lax-Oleinik formula 
cannot be used, and the correct weak solution can have shocks prohibited by entropy 
criteria (Section 3.1). Without differentiability of / we loose some good properties 
of characteristics (Section 3.2). In section 3.3 we prove Theorem 2, the uniqueness 
result. 

Section 4 reviews the essentials of the variational coupling approach to exclusion 
and zero-range processes from Seppalainen (1998a,b,1999), and develops the for- 
mula for the second class particle. Section 5 proves the laws of large numbers for 
the second class particle. The proofs are an interplay of the particle level properties 
from Section 4 and the properties of the macroscopic characteristics. 

Notation. The set of natural numbers is N = {1, 2, 3, . . . }, and the set of non- 



2. Results 

2.1. Hydrodynamic limits and entropy solutions. We start by reviewing the 
hydro dynamic hmit of totally asymmetric iiT-exclusion. The parameter i^ is a fixed 
positive finite or infinite integer iiT G N U {oo}. Each process consists of indistin- 
guishable particles that move on the one-dimensional integer lattice Z, subject to 
the restriction that a site i G Z can be occupied by at most K particles. If i^ = oo 
this constraint is absent. The state space of the process is iS = {0, 1, 2, . . . , K}"^ 
for finite X, and S = Z^ if K = oo. A generic element of S is rj = {rj{i) : i G Z), 
where r]{i) denotes the number of particles present at site i. 

The dynamics is such that at exponential rate 1, one particle moves from site i 
to z + 1, provided there is at least one particle at i and at most K — 1 particles at 
i + 1. The jump attempts happen independently at all sites i. 

To realize this dynamics, let {Di : z G Z} denote a collection of mutually inde- 
pendent rate 1 Poisson point processes on the time axis (0, oo). We usually think of 
-Dj as the random set of jump times (or epochs, to borrow Feller's term) rather than 
as the corresponding random step function. At each epoch of Di a jump from site 
i to site i -|- 1 is attempted, and successfully executed if the rules permit. The state 
of the process at time t is the configuration r]{t) — {r]{i, t) : i E Z) of occupation 
numbers. 

The dynamics can be represented by the generator L that acts on bounded cylin- 
der functions / on S: 

(2.1) Lfir^) = J2 Hvi^) >1,V{^ + 1)<K- l}[f{v'''+') - fiv)]- 

Here ?y*'*+^ is the configuration that results from the jump of a single particle from 
site i to site i + 1: 

r]{i) - 1, j = i 

V'''^Hj)=\ r?(^ + !) + !, i = ^ + l 

Special cases, (a) When K = 1 this process is known as the totally asymmetric 
simple exclusion process (TASEP). The Bernoulli probability measures {i'^} on S, 
indexed by density p G [0, 1], are invariant for the process. These are defined by 
the condition 

^P{ri{t^) = ^{i^) = ... = n{i^) = 1} = p-^ 

for any finite set {zi, . . . , ij^} of sites. 

(b) If i^ = oo this process is a member of the family of totally asymmetric 
zero-range processes (TAZRP). The i.i.d. geometric probability measures {i'^}, p G 
[0, oo), are invariant. These are given by 

(2.2) i^'iriih) = h,7]ii2) = A:2, . . . , vi^m) = k^} = H j^ (l^j 

for distinct sites ii, . . . ,im and arbitrary integers /ci, . . . , /c^n ^ Z_|_. 

(c) For 2 < i^ < oo we call this the totally asymmetric X-exclusion process. For 
this case the existence of spatially ergodic invariant measures for all densities has 



The basic result about the large-scale behavior of these processes is the following 
hydrodynamic scaling limit. Let po be a bounded measurable function on R such 
that < po{x) < K. Assume given a probability space on which are defined a 
sequence of initial configurations {r]n{i,0) : i G Z), n G N, and the Poisson jump 
time processes {Di}. The processes {?7n(')} ^^^ constructed on this probability 
space by the familiar graphical method. Assume that at time zero the following 
strong law of large numbers is valid: 



(2.3) For all — oo < a < 6 < oo, lim — N ?7^(z,0)= / pQ{x)dx a.s. 



Theorem 1. Under assumption (2.3) a strong law of large numbers holds for all 
t > and all — oo < a < 6 < oo.- 



(2.4) hm — y^ r]n(i,nt) = / p(x,t)dx a.s. 



i=[na] + l 

The deterministic limiting density p{x,t) is a weak solution of the Cauchy problem 

(2.5) Pt + /(p)x = 0, p{x,Q)=po{x). 

The flux function f of the conservation law (2.5) is a fixed, deterministic concave 
function on [0, i^] that depends on K. 

The limit (2.3) is the only assumption needed in the following strong sense: The 
joint distribution of the initially defined random variables {(?7n(^, 0) : i G Z) : n G 
N}, and {Di : z G Z} may be arbitrary. As long as the marginal distribution of 
r]n{i, 0) satisfies (2.3), and marginally {-Di} are i.i.d. rate 1 Poisson point processes, 
Theorem 1 is valid. 

A second point is that if (2.3) holds in probability, then the conclusion (2.5) also 
holds in probability. 

For exclusion [K = 1) and zero-range {K = oo) Theorem 1 has gone through 
many stages of generalization, beginning with the seminal work of Rost (1981). 
The entropy technique of Rezakhanlou (1991) proves it for multidimensional exclu- 
sion and zero-range processes, without requiring total asymmetry. The variational 
coupling method of Seppalainen (1998a,b) permits the most general initial distribu- 
tions for the one-dimensional, totally asymmetric case. For 2 < K < oo Theorem 
1 was first proved in Seppalainen (1999). 

Let us write fx for the fiux function / in (2.5) when we emphasize dependence 
on K. fx is known in cases where invariant distributions are explicitly known: for 
exclusion /i(p) = p{l — p) for < p < 1, and for zero-range /oo(p) = p/(l + p) 
for < p < oo. For 1 < i^ < oo we know the concavity and continuity of /k, 
the symmetry fxip) = fK{K — p), the monotonicity fxip) < /i^+i(p), and these 
bounds: 

(2.6) min|p(l - p) , 1/4} < /^(p) < — ^ for < p < K/2. 



Currently we have no information about the regularity of / for 1 < i^ < oo. 
As suggested by /i and /oo, we would expect / to be C^ and strictly concave, 
in other words, for its graph to have no corners or line segments. This lack of 
information about / creates a problem because the common definitions of entropy 
solutions of (2.5) are valid only for strictly convex or concave fiux functions. Thus 
we cannot apply standard entropy criteria such as one-sided bounds on jumps or 
the Lax-Oleinik formula. We show this in Section 3. 

The next task is to complement Theorem 1 with a uniqueness criterion. First we 
explain how the solution p{x,t) that the particle system chooses in Theorem 1 is 
defined. In the next discussion assume first K < oo. Afterwards we comment on 
the K = oo case. 

Define an antiderivative uq of po by 

(2.7) uo{0) = and uo{b) — uo{a) = / po{x)dx for all a < b. 



Let g be the nonincreasing, nonnegative convex function on R defined by 

(2.8) g{x)= sup {f{p)-xp}. 

0<P<K 

Symmetry of / and bounds (2.6) imply that /'(O) = 1 and f'{K) = —1, the 
existence of the derivatives being part of the conclusion. It follows that g satisfies 

(2.9) g{x) = for a; > 1 and g{x) = —Kx for x < —1. 
Then for a; G R set u{x^ 0) = uq{x), and for t > 

' X — 



(2.10) u{x,t) = snY>\uQ{y)-tg 

yen L 

The supremum is attained at some y E [x — t, x + t]. The function u is uniformly 
Lipschitz on R X [0, oo), nonincreasing in t and nondecreasing in x. (2.10) is known 
as the Hopf-Lax formula. It defines u{x, t) as a viscosity solution of the Hamilton- 
Jacobi equation 

(2.11) wt + /(w,) = 0, u{x,Q)=uo{x) 

[Ch. 10 in Evans (1998)]. The viscosity solution is unique by Theorem 2.1 of Ishii 
(1984) that applies to (2.11) because / is continuous and u is uniformly continuous. 

If i^ = oo we know that g{x) = (1 — \/xY for < x < 1, (7 = on (1, 00) 
and we can take g = 00 on (— oo,0). The supremum in (2.10) is attained at some 
y G [x — t,x\. Since po is assumed bounded, u is again Lipschitz on R x [0, 00). 

Let 

(2.12) p{x,t) = —u{x,t) 

be the a.e. defined x-derivative. Recall this standard definition: a locally bounded 
measurable function A(x, t) is a weak solution of (2.5) if for all compactly supported 
and continuously differentiable test functions (/> on R x [0, 00), 



(2.13) / / \\(xA)(j)t{x,t)-\-f{\{xA))(t)Ax,t)]dxdt-\- / pn{x)(j){x,Q)dx 



0. 



By the proof of Thm. 3.4.2 in Evans (1998), p{x,t) is a weak solution of (2.5) in 
the integral sense (2.13). The solution p{x,t) defined by (2.12) is the one chosen 
by the particle system in Theorem 1. 

In cases K = 1 and K = oo we know that the flux / is differentiable and strictly 
concave. Then p{x, t) can be equivalently defined by the Lax-Oleinik formula, and 
we can call p{x,t) the "unique entropy solution." For 2 < i^ < oo we can only 
define p{x,t) through (2.12). In Sections 3.1-3.2 we show that without further 
knowledge about / this solution cannot be defined by the Lax-Oleinik formula, and 
we cannot rule out discontinuities of the wrong kind that violate entropy conditions. 

To complement the characterization of p{x,t) through the viscosity solution of 
(2.11), we give a uniqueness criterion for which strict concavity or differentiability 
of / are not needed. 

Theorem 2. Among the weak solutions o/ (2.5), p{x,t) o/ (2.12) is characterized 
by maximal current over time. More precisely, suppose A(x, t) is a nonnegative, 
locally hounded measurable function on Rx (0, oo) that satisfies the integral criterion 
(2.13). Fix t > 0. Then for almost all x eH, 

(2.14) / f{X{x,s))ds< [ f{p{x,s))ds. 

Jo Jo 

Conversely, if equality holds in (2.14) for almost all (x^t), then A(x,t) = p{x^t) for 
almost all (x, t). 

Thus the particle system chooses the weak solution that transports the most 
material. Theorem 2 is proved in Section 3.3. The system maximizes the current 
because / is concave. With a convex flux the current would be minimized, as is the 
case with the stick model version of Hammersley's process, see Seppalainen (1996). 

2.2. Second class particles and characteristics. In this section we state the 
laws of large numbers for the second class particle. The representation of the second 
class particle in the variational coupling is in Section 4.2. 

Let X{t) denote the position of a single second class particle added to the process 
r]{t). X{t) is rigorously defined as follows. Assume given a probability space on 
which are defined the following random variables: the Poisson jump time processes 
{Di}, an initial particle configuration {ri{i, 0) : z G Z), and an initial location ^(0). 
Define a second initial configuration (77(2, 0) : i G Z) that differs from rj only at site 
X(0): ?7(X(0), 0) = ?7(X(0), 0) + 1 and ?)(«, 0) = r]{i, 0) for i ^ X(0). [Of course we 
must have r]{X{Q), Q) < K — 1 a.s. for this to be possible.] 

Run the processes r]{t) and f}{t) so that they read the same Poisson jump time 
processes {Di}. This is known as the basic coupling of ?7(t) and fi{t). There is always 
a unique site X{t) at which the two processes differ: fi{X{t),t) = r]{X{t), t) + l and 
fi{i,t) = rj{ijt) for i 7^ X{t). This defines the position X{t) of the second class 
particle. 

The way in which the second class particle yields to first class particles becomes 
evident when we consider the two types of transitions that can happen to it. Let 

X{T-)=i. 

(i) Suppose r is a jump epoch for Poisson process Di. The second class particle 
can jump from i to i + 1 if there is space at site i + 1, but only if there is no first class 
particle to jump. In other words, if rj{i, r— ) = and r]{i + l, r— ) < K, then X{t) = 



(ii) Suppose r is a jump epoch for -Di_i. The second class particle jumps from 
i to z — 1 if first class particles fill site i. In other words, if r]{i — 1, r— ) > 1 and 
?7(z, T-) = K - 1, then X{t) = i - 1. 

If K = oo the second class particle never jumps left and only the first kind of 
transition happens to it. 

The definition in terms of ry and fj generalizes naturally to a two-species process 
of first and second class particles. But in this paper we consider only a single second 
class particle added to the process. 

Next we discuss the characteristics of the conservation law (2.5). The lack of 
information about the regularity of / becomes a problem again. The well-known 
description of characteristics as generalized solutions of the ordinary differential 
equation dx/dt = f'{p{x,t)) [see Dafermos (1977)] does not work now. As Reza- 
khanlou (1995), we use the Hopf-Lax formula to define characteristics. 

Once u{x, t) is defined by (2.10), a semigroup property is in force: for all < s < t 
and a; G R, 

I I X — V 

(2.15) u{x, t) = sup <^ u{y, s) - {t - s)g 

yen[ \t-s 

The supremum is assumed at some point y E [x — {t ~ s) , x + {t — s)]. (If K = oo 
this range would actually be [x — {t — s), x].) By continuity the supremum in (2.15) 
is attained at the least and largest maximizers, denoted by 

I f X — V 

(2.16) y~{x; s,t) = inf < y > X — {t — s) : u{x, t) = u{y, s) — {t — s)g 



J-s 
and 

■^ (T-Mf.) ^ fiun < 11 <C T -\- (f — .si • i/Y.r, f.) = iidr fi\ — (f. — fi\n\ 

t- S 



(2.17) y {x; s, t) = sup \y<x + {t — s): u{x, t) = u{y, s) — {t — s)g 



Abbreviate y^(x, t) = y^{x; 0, t). As functions of x for fixed < s < t, y~^{x; s, t) is 
right-continuous, y~{x; s,t) is left-continuous, both are nondecreasing, and trivially 
from the definition lim^j^-oo y^(a^; s, t) = — oo and lim^j^oo y^(3^; s, t) = oo. See 
comment 2.3.1 below for justification of the restriction x — {t — s) < y < x + {t — s) 
in the definitions of y^(x; s, t). 

The minimal and maximal forward characteristics are defined for 6 G R, < s < 
6, as 



(2.18) x~{b;s,t) =mf{x ■.y+{x;s,t) > b} 
and 

(2.19) x^{b; s, t) = sup{a: : y~{x; s, t) < b}. 
Immediate properties are these: 



10 

As functions of 6, x^{b;s,t) is right-continuous, x~{b]Sjt) is left-continuous, and 
both are nondecreasing. For the s = case abbreviate a;^(6, t) = x^(6;0, t). As 
functions of t, x^(6; s, t) are Lipschitz and satisfy the initial condition x^(6; s, s) = 
b. When characteristics are unique we write x{b] s, t) = x^{b; s, t). 

Now consider the setting of the hydrodynamic limit Theorem 1, with a sequence 
of processes rjn that satisfy the initial limit (2.3). Let Xn{t) denote the position 
of a second class particle in the process r]n- The initial position Xn{0) may be 
deterministic or random, and may depend in an arbitrary fashion on the initial 
particle configurations {r]n{i,0)) and the Poisson processes {Di}. For the first 
result assume a law of large numbers at time zero: for a deterministic point 6 G R 

(2.20) lim ^^ = b a.s. 

Theorem 3. Under assumptions (2.3) and (2.20), for any fixed t > 0, 

(2.21) X- {b, t) < lim inf ^^^ < lim sup ^^^ < x+ {b, t) 

almost surely. In particular, if x{b^ t) = x^(6, t), then we have a strong law of large 
numbers: 

(2.22) lim ^IiM=^(6,t) a.s. 

In the special case where characteristics are unique we can dispense with assump- 
tion (2.20). The second class particle tracks the macroscopically defined charac- 
teristic that emanates from its random initial position, as stated in part (a) of the 
next theorem. We only need to assume that the fiuctuations of the initial location 
are macroscopically bounded as n grows: 

(2.23) — cxD < lim inf — - — < lim sup — < oo almost surely. 

The uniqueness of characteristics is true if we look at characteristics started at a 
time s > 0. But only if the fiux function is differentiable, so this assumption needs 
to be made explicitly for 2 < K < oo in part (b). 

Theorem 4. Assume (2.3) and (2.23). 

(a) If x{y,t) = x^{y,t) for all y eH, then almost surely 



(2.24) hm 



>CXD 



X^jnt) ^ r xr^jo) ^ 



n \ n 



0. 



(b) Let < s < t. In the case 2 < K < oo suppose f is differentiable on [0, K]. 
For K = 1 and K = oo this is known and does not need to be assumed. Then 
x{y; s, t) = x'^{y; s, t) for all y E H, and almost surely 



(2.25) lim 



*oo 






n \ n 



0. 



The meaning of x (n~^X^(0) , t) is that the deterministic function x{-,t) defined 
by (2.18-19) is evaluated at the random point n~^X^(0). If assumption (2.23) is 
strengthened to na < X„(0) < nb for all large enough n, then for Theorem 4(a) it 
is sufficient to assume x{y,t) = x^{y,t) for a < y < b. As for Theorem 1, if the 
assumptions are valid in probability, then so are the conclusions. 
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2.3. Comments on the assumptions and results. 

2.3.1. Differentiability and strict convexity. Since g is convex, left and right 
derivatives 

// X ,. 9{x-e)-g{x) '( \ v 9{x + e) - g{x) 
a (x) = Iim and Q^(x) = lim 

exist at all points x. Same is true for /, and for both / and g differentiability on an 
open interval is equivalent to continuous differentiability. Strict convexity means 
that there are no xi < X2 such that g'_^_{xi) = g'_{x2), which is the same as not 
having a line segment in the graph of g. From the duality (2.8) it follows that / is 
differentiable on [0, K] iff g is strictly convex on [—1, 1], and / is strictly concave on 
[0, K] iff g is differentiable on R. Consult Rockafellar (1970) for a general treatment 
of convex analysis. 

Theorem 4(b) is true because strict convexity of g guarantees that x{y; s, t) = 
x^{y; s,t) for all y G R, if s > 0. This is proved in Proposition 3.1 in Section 3.2. 

The restriction x — {t — s) < y < x + {t — s) in definition (2.16-17) of y^(x; s, t) 
is necessary because g is linear outside [—1, 1]. For example, take as initial profile 
Po{x) = K for X < and po{x) = for x > 0. Then for x > t, the entire interval 
[0, X — t] consists of maximizers in the Hopf-Lax formula. Without the restriction 
in (2.16) we would have y~{x^ t) = for all x > t, and then a;+(0, t) = oo which is 
not the appropriate characteristic. 

2.3.2. Shocks. Consider the simplest shock case (Riemann problem) with initial 
profile po{x) — al{x < 0} + (31{x > 0} where a < (3. The solution (2.12) is 
p{x,t) = al{x <^t}+(31{x > ^t} with shock speed ^ = {f{(3)-f{a))/{(3-a). The 
unique characteristic from the origin is a;(0, t) = ^t. Theorem 3 says that the second 
class particle converges to the shock. This generalizes results of Ferrari (1992) in 
two ways: From exclusion to X-exclusion, and from initial product distributions to 
general initial distributions that have the shock profile. 

2.3.3. Rarefaction fan. Ferrari and Kipnis (1995) proved convergence in prob- 
ability in (2.25) for the following case of TASEP: initial distribution is product 
measure with density p to the left and density A to the right of the origin, with 
p > A, and the second class particle initially at the origin. This initial pro- 
file po{x) = pl{x < 0} -f- Al{a; > 0} has an inadmissible shock at the origin, 
which produces a "rarefaction fan" of characteristics with x~{0, t) = (1 — 2p)t and 
x~^{0,t) = (1 — 2A)t. Ferrari and Kipnis also proved that weakly n~^Xn{nt) con- 
verges to the uniform distribution on [x~(0, t), x+(0, t)]. Consequently one cannot 
hope for a law of large numbers when x~{0, t) < x^(0, t). 

If / is not differentiable at p, then the situation x~{b,t) < x^{b,t) happens even 
for the constant profile p{x, t) = p. See Example 3.1 in Section 3.2. In this case we 
cannot deduce convergence of n~^Xn{nt) from Theorem 3, only the bounds (2.24), 
unless we can prove the differentiability of /. 

2.3.4. The characteristic ordinary differential equation. When / is strictly con- 
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x{t) — a;^(6, t) are Filippov solutions of the o.d.e. 

(2.26) - = f'{p{x,t)),x{0)=b, 

where p{x,t) is the Lax-Oleinik solution of (2.5). Since p{x,t) does not exist at all 
X, (2.26) is interpreted as 

dx 

(2.27) essinf /'(p(x,t)) < — < esssup /'(p(x, t)) for a.e. t. 

^ at X 

The essential infimum relative to x, at the point (x, t), is defined by 
essinf/'(p(x,t)) = lim essinf f'{p{y,t)), 

X 5\0 x—S<y<x+S 

and similarly for the essential supremum. The Lax-Oleinik formula (we review it 
in Section 3.1 below) shows that (2.27) is equivalent to 

(2.28) /' (-g' f^^yli^]] <^<f' f-g' ("^ - ^"(^'^) 



t J J - dt - \ \ t 

for almost every t. These results are developed in Rezakhanlou (1995) for the 
exclusion and zero-range processes. 

At present we cannot ascertain strict concavity or differentiability of / for 2 < 
K < oo. This is a problem, because to interpret (2.26) in the Filippov sense, 
f'{p{x,t)) must exist at least a.e. But if /' does not exist at p and we consider the 
constant solution p{x,t) = p, then f'(p{x,t)) does not exist at any (x,t). 

Without any assumptions on / beyond concavity, bounds (2.28) are valid for 
x{t) — x^{b, t) in the form 

(2.29) /; (-si (i^i^) ) < I < /i [-,'. ( "-";'"■" 

for almost every t. If we assume / differentiable, then in fact (2.26) is valid again 
in the sense (2.27), with p{x,t) given by (2.12). In other words, the o.d.e. works 
even though the Lax-Oleinik formula fails without strict concavity of /. We omit 
the proof of (2.29) and the last statement above. 

2.4. Generalization: particles jump in batches. We discuss here a generaliza- 
tion of the i^-exclusion and zero-range processes for which the variational coupling 
technique works and Theorems 1-4 are true. 

From a queueing perspective, the particles are customers moving through an 
infinite sequence of servers. Each server has space for K customers in its queue, and 
r]{i,t) is the number of customers waiting at server i at time t. In the formulation 
(2.1), each server serves at exponential rate 1, and only one customer is served at a 
time and moves on to the next server. We generalize this by permitting customers 
to move in batches of varying sizes, with exponential rates and independently for 
each batch size. Let f3h be the rate at which a batch of h customers can jump from 
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customers to jump [so either r]{i) < h or rj{i + 1) > K — h], we do not suppress the 
entire jump but instead move as many customers as possible without violating the 
constraints. 

The family of processes is parametrized by K and a sequence of nonnegative 
numbers {(3h : /i G N). If i^ < oo assume (3h = for h > K. To realize this 
dynamics take a collection of mutually independent Poisson point processes {Df : 
i E Z,h E N} on the time axis (0, oo), so that the rate of D^ is (3h- At each epoch 
of -Df , a batch of h particles (customers) attempt to jump from site i to site i + 1. 
The generator becomes 

where t^'^'^-'+i is the configuration that results from an attempt to move a batch of 
h particles from site i to i + 1. Denote the actual number of particles that move by 

b = b{hj r]{i), r]{i + 1)) = minj/i, ?7(z), K — r]{i + 1)}. 

Then the new configuration is 



h,i,i+l/j\ 



The definition of a second class particle is the same as before. Prom the queueing 
perspective it is a second class customer. He can move forward only if first class 
customers do not fill the batch that is served, and he has to move back to the 
previous server if arriving first class customers fill the capacity of the queue. 

The hydrodynamic limit Theorem 1 is valid as stated for the batch process too, 
provided the jump rates are restricted so that we can expect finite limits. It is 
enough to require X^/ieN/^^^^ ^ ^^- Under this assumption a weaker version of 
Theorem 1 with convergence in probability was proved in Seppalainen (2000) . That 
paper also shows how to adapt the variational coupling approach to batch jumps. 

li K = oo, the i.i.d. geometric probability measures (2.2) are again invariant, for 
any choice of the rates (/3/i). In this case the fiux function is 



r,{{) - b, 


j =i 




T]{i + l)+b, 


j =i 


+ 1 


vij), 


3+^ 


,i + l. 



/(p) = E ^^ 



/i6N 



' +^r^ U---+^ ' 



1 + p Vi + py Vi + p 



3. Scalar conservation laws with concave flux functions 

In this section we investigate the consequences of the lack of strict concavity and 
differentiability of the fiux function. It is important to clarify these facts because 
they directly infiuence the results we obtain for the particle models. Section 3.1 
looks at strict concavity and section 3.2 at differentiability. In section 3.3 we prove 
the uniqueness criterion of Theorem 2 which assumes neither strict concavity nor 
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3.1. Strict concavity of the flux and entropy shocks. We begin by recalling 
the Lax-Oleinik formula from Section 3.4 in Evans (1998), and then show how it 
depends on the strict concavity of the flux function. 

Lax-Oleinik formula. Suppose / is strictly concave and differentiable on its domain 
[0,iir]. Deflne two functions 

(3.1) p^{x,t) = -g \ ^ 

For flxed t > 0, under these assumptions y^{x,t) = y~{x,t) for all but countably 
many x, and consequently the Lax-Oleinik solution 

(3.2) pi^o{x,t) = p^{x,t) 

is well-deflned for all but countably many x. It satisfles p-Lo{x,t) = Ux{x,t) a.e., 
and is a weak solution of (2.5) in the integral sense (2.13). 

In particular, when it exists plo is the solution deflned in (2.12) for the hydro- 
dynamic limit. One can check that Plo(' ^t) is continuous at all but the countably 
many x at which p~{x,t) ^ /9+(x,t). If /Olo('5^) is discontinuous at x, there must 
be a jump upward: p-Lo{x—,t) < p-Lo{x+,t). This is because g' is nondecreasing 
and y~{x, t) < y^{xj t). Let us call such a jump an entropy shock for a concave flux 
function. [For a convex flux entropy shocks are downward jumps.] 

There are two ways deflnition (3.1)-(3.2) could conceivably run into a problem: 

(i) If / is not strictly concave then g' fails to exist at some points and the right- 
hand side of (3.1) may not be deflned. 

(ii) If /' does not exist everywhere then y~{x,t) < y^{x,t) may happen for a 
nontrivial interval of points x, and formula (3.2) is not true a.e. See Example 3.1 
in Section 3.2. 

It turns out that if/ is strictly concave, deflnition (3.1-2) works even if y~(x, t) = 
y^{x,t) fails. And all discontinuities are entropy shocks. 

Theorem 5. Suppose the flux f is strictly concave. Then for a fixed t > the 
solution p{x,t) = Ux{x,t) defined as in (2.12) exists and is continuous on a set 
H C H that contains all but countably many x. On H p{x,t) — p^{x^t) so the 
Lax-Oleinik formula is valid. For all x p{x—,t) < p{x+,t), so the discontinuities 
of pi' J t) o,re entropy shocks. 

Conversely, suppose f has a linear segment. Then there is a range of densities p 
such that the constant solution p{x,t) = p cannot be represented by any left or right 
derivative of g, so the Lax-Oleinik formula (3.1-2) fails. Furthermore, the relevant 
solution defined by (2.12) can have nonentropy shocks p{x—,t) > p{x+,t). 

Proof. Strict concavity of / implies that g is differentiable. Inequalities (3.3)-(3.4) 
below then imply that the right and left ^-derivatives of u{x, t) exist and are given 

by 

/ ,N I ( X — y {x,t)\ ^ 
Ux±[x,t) = -g I ^ \=p [x,t). 

Since g' is continuous and y^{- , t) nondecreasing and right-continuous, we see that 
/9"'"(-,t) is right-continuous, has at most countably many discontinuities, and has 
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By the Lipschitz property w(x, t) is differentiable at a.e. x and is the integral of 
its derivative. Thus 

X fX 

±1 



Jo Jo 

Consequently p{x,t) = Ux{x,t) exists and equals p^{x,t) at any x where one 
of p^{x^t) is continuous. This defines the set H. From this follows p{x—,t) — 
p~{x,t) < p+{x,t) = p{x+,t). 

For the converse, suppose f'{p) = ^ for po < P < Pi- Then g has a corner at ^, 
with g'_{0 = ~Pi ^^d g'^iC) = —po- The values p G {po. Pi) are not taken by g'^ 
at any point, so a constant solution p{x,t) = p G (po^Pi) cannot be represented by 
a derivative of g. 

The initial profile ^0(2;) = K ioY x <Q and po{x) = for x > gives the solution 
u{x,t) = —tg{x/t). Then p((t^) — ,t) = pi > po = ^((^0+?'^) is a nonentropy 
shock, n 

To finish this section we prove the inequalities used in the proof above. Let us 
write lim to simultaneously include both the limsup and the liminf as £ \ 0. 

Lemma 3.1. The following inequalities hold for right and left x-derivatives of 
u{x, t) at fixed t > 0. 

, fx-y-{x,t)\ -^u{x-e,t)-u{x,t) , fx-y-{x,t) 

(3.3) -g^ I ) < Inn < -g. 

and 



e\o ^^ \ t 



, fx-y+{x,t)\ -^u{x + e,t)-u{x,t) , fx-y+{x,t) 

(3.4) -g^ < hm < -g 



^ J e\0 ^ \ t 

Proof. The proofs use convexity, the Hopf-Lax formula, and the right [left] conti- 
nuity of y+(- , t) and g'_^_ [y~{- , t) and g'_]. As an illustration, here is the argument 
for the second inequality of (3.4). We leave the rest to the reader. 

u{x + £, t) — w(x, t) 

, I , ,, f X + e — y^(x + e,t)\ 
< uo{y+{x + £, t)) - tg ( Y J 

-uo{y+{x + e,t))+tgi '^ ^j 

^_£./^iV/ f x-y^{x + e,t) e\ / x - y+{x + e,t) \\ 

, f x-y+{x + e,t) \ 

- ^-^H t )■ 

Consider two possibilities: If y^{x + £, t) = y~^{xj t) for small enough £ > 0, then 
, / x-y+{x + e,t) \ ^ , / x~y+{x,t) \ ^ , / x-y+{x,t) \ 
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The other possibihty is that y~^{x + £,t) > y~^{x,t) for all £ > 0. From right 
continuity of x \-^ y^{x,t) it follows that {x — y^{x + e,t))/t increases strictly to 
{x — y^{x, t))/t as £ \ 0. And then 

, [x-y+{x + e,t)\ , fx-y+{x,t) 
lim g, = fi*- 

In either case we obtain the second inequality of (3.4). D 

3.2. Differentiability of the flux and characteristics. In Corollary 3.1(b) 
below we show that differentiability of / implies y~{x,t) = y~^{x,t) for all but 
countably many x, the property used in part (3.2) of the Lax-Oleinik formula. 
In Proposition 3.1 we prove the property needed for Theorem 4(b), namely that 
x~{b;s,t) = x~^{b;Sjt) for < s < t. Example 3.1 shows how various properties 
can fail without differentiability of /. The hypothesis of differentiability of / is 
used in the equivalent form that g is strictly convex on [—1, 1]. 

Let I{x; s, t) denote the set of y G [x — {t — s),x + {t — s)] at which the supremum 
in (2.15) is achieved. 

Lemma 3.2. Suppose < ti < t2 < ts, y2 G I{x;t2,ts), and yi G /(2/2;^i, ^2)- 
Then yi G I{x;tijts). And either 

x-yi ^ x-y2 ^ y2-yi 

ts — tl ts — t2 t2 — ti 

or g is linear between (1/2 — yi)/(^2 — ^1) and {x — y2)/{ts — ^2)- 

Proof. Modify the proof of Rezakhanlou's (1995) Lemma 3.1. D 

Lemma 3.3. Suppose g is strictly convex on [—1,1], xi < X2, yi G I{xi;s,t) for 
i = 1,2. Then yi < 2/2- 

Proof. The assumptions yi G I{xi;s,t) give 

u{yi, s)-{t- s)g i ^_^^ j > u{y2, s) - {t - s)g i ^_^^ 

and 

u{y2, s)-{t- s)g ( -j—^ 1 > u{yi, s) - {t - s)g i ^_^^ 

To get a contradiction, suppose yi — y2 > 0. Then the above inequalities give 



< 

X2-y2 \ _ ( ^2-yi \ ( xi-y2 \ _ ( xi-yi 

t-s J \ t-s J y t-s J \ t-s 

Since X2 - y2 > xi - 1/2, X2 - yi > xi - yi, and {xi - yi)/{t - s) e [-1, 1], this 
contradicts the assumption on g. D 



17 

Corollary 3.1. Suppose g is strictly convex on [—1, 1]. 

(a) For xi < X2, y^{xi] s,t) < y~{x2] s,t). 

(b) For fixed < s <t, y~{x;s,t) = y~^{x;s,t) for all but countably many x. 

(c) x~^{b; s, t) = mf{x : y~^{x; s, t) > b} and x~{b; s, t) = sup{x : y~{x] s, t) < b}. 

(d) If x~{b;s,t) < z < x+(6; s, t), then y^^^z; s,t) =b. 

Proof, (a) is immediate from Lemma 3.3. 

(b) Suppose xq is a continuity point of the nondecreasing, right-continuous func- 
tion X i-H^ y^{x] s, t). Then 

y+(xo;s,t) > 2/"(xo;s,t) > hm y+(x;s,t) = y+(xo;s,t), 

where the second inequahty follows from (a). Thus y~{x;Sjt) = y~^{x;s,t) at all 
continuity points of y~^{x; s,t), which includes all but countably many x. 

(c) Let w{b) = mf{x : y~^{x;Sjt) > b}. We have x~^{b;s,t) > w{b) without any 
new assumptions: if a: < w{b), then 

y~^{x;Sjt)<b =^ y~{x;s,t)<b =^ x^{b]Sjt)>x. 

For the converse we need part (a). Let x > w{b). Then there exists xi such that 
X > xi > w{b). 

y^{xi;s,t)>b =^ y~{x]s,t)>b =^ x^{b;s,t)<x. 

We leave the proof of the other formula to the reader. 

(d) Immediate from definition (2.18-19) and part (c). D 

Proposition 3.1. Suppose g is strictly convex on [—1,1], and s > 0. Then 
x~{b] s, t) = x~^{b; s, t) for all b eH and t > s. 

Proof. Suppose there exists a z such that x~{b; s,t) < z < x~^{b; s, t). By Corollary 
(d) y^{z; s, t) = 6, and consequently for e > 0, 

u{b, s)-{t- s)g (^j > Hb ±e,s)-{t- s)g r ~_^J^ 

Letting e \0 then gives 

u(b + e,s) - u(b,s) ^ , fz-b 

hm sup < —g_ 

£\o £ \t - s 



and 



,. . u{b-e,s)-u{b,s) , fz-b 

limmt > —gj^ 



e\0 —e '^ \t — s 



Since we may vary z in an interval, and since by assumption the slope of g cannot 
remain constant in this interval, we conclude the strict inequality 

u(b + e,s) — uib, s) , ^ uib — e,s) — uib, s) 
limsup ^ —^ -^-- < hmmf ^ —^ -—--. 
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By (3.3-4), this implies 

(3.5) g^ ]> g. 



The slope of a convex function is nondecreasing and by definition {b~y~{b, s))/s > 
{b — y~^{b, s))/s. So (3.5) implies that y~{b,s) = y+(6, s) and g has a corner at 
(6 — y^{b^s))/s. Since we are assuming g strictly convex in [—1,1], Lemma 3.2 
now forces {z — b)/{t — s) = (6 — y^{b,s))/ s. This cannot hold for distinct z''s. 
But the preceding argument is valid for all z G (x~(6; s, t), x^{b; s, t)), so to avoid 
contradiction we must have x~{b] s, t) = x^{b; s, t). D 

Finally we show by example how a point of nondifferentiability in / destroys all 
the good properties of the characteristics. 

Example 3.1. Suppose p G (0, K) is a point such that f'_{p) = ^2 > ^i = f+{p)- 
Then g has a linear segment of slope — p over [^1,^2]- Apply the formulas to the 
constant initial profile po{x) = p, with ^0(3^) = px- We get I{x;s,t) = [x — 
^2{t - s),x - ^i(t - s)], y~{x]s,t) = X - ^2{t - s), y^{x]s,t) = x - ^i(t - s), 
x~{b; s, t) = 6 + ^i(t — s), x~^{b; s, t) = b + ^2{i — s), and u{x, t) = px — tf{p). 
In particular p{x^t) = Uj:{x,t) = p, so the constant profile is the solution. But 
the conclusions of Lemma 3.3, Corollary 3.1(a)-(b), and Proposition 3.1 have been 
contradicted. 

3.3. Proof of Theorem 2. We prove the uniqueness criterion of Theorem 2 under 
these general assumptions: The fiux / is continuous and concave on an interval /, 
and g is defined by g{x) = suppgj{/(p) — px}. The initial density po{x) is a given 
locally bounded measurable function on R, u{x, t) is a locally Lipschitz function on 
R X [0, 00) defined by the Hopf-Lax formula (2.10), and p(x, t) = u^i^x, t) exists for 
a.e. X for any fixed t. A(x, t) is a locally bounded measurable function on Rx [0, 00) 
that satisfies the boundary condition A(a;, 0) = po{x) a.e. and the integral condition 
(2.13). The ranges of p and A lie in /. 
The goal is to show that, for a.e. x, 

(3.6) / f{X{x,s))ds< [ f{p{x,s))ds 

Jo Jo 

for all t > 0. And conversely, if equality holds in (3.6) for almost all {x,t), then 
A(x,t) = p{x,t) for almost all (x,t). 

Some technical difficulties stem from assuming only (2.13) and no regularity on 
A. We cannot expect the function x ^^ A(a;, t) to be a sensible object for every fixed 
t because A(-, t) can be redefined on any null set of t's without affecting (2.13). At 
t = we have 

(3.7) lim - I dt I X{x,t)ijj{x)dx = / po{x)iJj{x)dx 
e^o £ Jo Jr Jr 

for i/j G C^(R). This is obtained from (2.13) by taking (/){xjt) = 'i/j{x)g{t) where 
(7(0) = 1, g' = —l/e on (0, e), and g{t) = for t > e. [The C^ test functions in 

('o 1o^ „. — 1 1„„ — ] u,. „ — i-: „ „; — .„„.;„„ r^l j-„„j- t, .j-; — „ u,. j-„u; — , ^:,^:^-„ 1 
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A corresponding property holds for certain later times t. Abbreviate 

Leb(/i) = the set of Lebesgue points of the function h 

for any measurable, locally bounded function h on either R or R x [0, oo). [For 
the definition and properties of Lebesgue points see Section 7.2 in Wheeden and 
Zygmund (1977).] Define 

T = {t > : (x, t) G Leb(A) for almost every x G R}. 

Since the complement of Leb(A) is a null set of R x (0, oo), almost every t lies in 
T. By considering test functions of the type (j){x,t) = 'i/j{x)g{t) for suitable g, one 
obtains from (2.13) a limit such as (3.7) around t G T, and then the following for 
t G T: 

(3.8) / iJj{x)X{x,t)dx — / ip{x)po{x)dx = / / f{X{x,s))ip'{x)dxds 

Jr Jr Jo Jr 

for compactly supported C^ functions V' on R. By taking limits of such functions, 
(3.8) is valid for compactly supported, continuous, piecewise C^ functions i/j. 

Lemma 3.4. (i) Fix finite A < B and T > 0. Then there exists a constant 
C = C{A,B,T) such that, for all A < a < b < B and s,t e T U {0} such that 
0<s,t<T, 



(3.9) 



b j-b 

\{x,t)dx — I \{x,s)dx 



< C ■\t-s\. 



(ii) // 1 G T U {0} and (a^t), {h,t) G Leb(/ o A), then we have a time derivative 
along points of T: 

(3.10) lim^— i / \{x,t)dx- I X{x,s)dx} = f{X{a,t)) - f{X{b,t)). 

Proof. In (3.8) take i/j such that i/j = outside [a — e, 6 + e], i/j' = 1/e on (a — £, a) 
i/j = 1 on [a, 6], and i/j' = —1/e on (6, b + e). We get 

,.& pb 

\{x,t)dx — I X{x, s)dx + 0{e) 

(3.11 



-1 pi pa -| pz i*o-\-s 

- / f{X{x,T))dxdT / / f{X{x,T))dxdr. 

£ Js Ja-e £ Js Jb 



The 0{e) term contains the integrals over [a — e, a] and [6, b + e] from the left-hand 
side. These are bounded by 0{e) uniformly over the range of a, 6, s, t by the local 
boundedness of A. The 1/e factors on the right are the slopes oft/j. The second line 
in (3.11) is bounded by C\t — s\ in absolute value, where C is twice the supremum 
of |/(A(x, r))| over [A, B] x [0, T]. Part (i) follows by letting e\0. For part (ii) set 

. xIj. „I j;,.;j„ j-U, ,u /'o i i ^ u,. j. „ ^„^- -C„„j- rn -> „ , j. . 1 i.i ;: \ n i — i 
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As a locally Lipschitz function u{x, t) is differentiable a.e. by Rademacher's the- 
orem [Section 5.8.3 in Evans (1998)]. And as the viscosity solution defined by the 
Hopf-Lax formula, u{xjt) satisfies the equation ut + f{ux) = a.e. [Theorem 5, 
Section 3.3.2 in Evans (1998)]. By definition p = Ux exists a.e., so for almost every 
xq and all t > 0, 

(3.12) / f{p{xo,s))ds = - ut{xo,s)ds = uo{xo) -u{xo,t). 

Jo Jo 

It is also the case that for almost every xq, (xq, t) G Leb(/ o A) for almost every 
t. Fix such an xq for which also (3.12) holds. For (a;,t) G R x T define 



(3.13) v{x,t)= / X{y,t)dy- / f{X{xo, s))ds + uo{xo). 

Jxo Jo 

The first integral is interpreted with a sign so that v{x2, t) —v{xi,t) = f ^ A(y, t)dy 
for xi < X2- By Lemma 3.4(i) and local boundedness of A and / o A, t; is locally 
Lipschitz on Rx T. Consequently v extends uniquely to a locally Lipschitz function 
on R X [0, oo). This extension has the correct boundary values at t = because by 
Lemma 3.4(i), 

pX pX 

lim v{x,t)= / X{y,0)dy + uo{xo) = / po{y)dy + uo{xq) 
^uo{x). 

The goal is now to prove 

(3.14) v{x,t) > u{x,t) for all (a;,t). 

Setting X = Xq in (3.14) then gives, by (3.12) and (3.13), the conclusion (3.6) for 
X = Xq and all t G T. Both sides of (3.6) are continuous in t because the integrands 
are locally bounded, so (3.6) follows for all t. This whole argument beginning with 
(3.12) can be repeated for almost every xq. At the end of this section we argue 
that equality in (3.6) implies X = p a.e. First we prove (3.14). 

Fix (x, t). For a fixed y G R set ^(s) = y + {s/t){x — y) for s G [0, t]. As y varies 
we consider the collection of line segments {(^(s),s) : < s < t} on the plane 
R X [0, oo). Their union forms the "backward light cone" from the apex (x,t) to 
the base line R x {0}. For almost every such y, it must be that 

{C{s),s) G Leb(/ o A) n Leb(A) 

for almost every s G [0, t], because only a null set of points of R x [0, oo) fail to be 
Lebesgue points. Fix now such a y. 

The function 7(s) = v{^{s), s), < s < t, is Lipschitz, and hence Y{s) exists for 
a.e. s and 

(3.15) 7(t)-7(0)= / 'y'(s)ds. 
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Fix So G [0,t] such that aU this holds: Sq G T, 7'(so) exists, (iCo, Sq)? (^(■§0)7 ■So) G 
Leb(/ o A) n Leb(A), and sq G Leb{/(A(xo, ■))}• These conditions are satisfied by 
a.e. So- For such an sq we prove that 

(3.16) Yiso) = ^A(e(so), so) - /(A(^(so), so)). 

To prove (3.16) we calculate 7'(so) as the limit 

s\so S — So 

SET 

So consider s G T, s > Sq. 

7(g) -7(^0) ^ v{C{s),s)-v{C{so),so) 

s - So s - So 

= / Hy,s)dy-\ <^ / X{y,s)dy- / X{y,so)dy} 

1 r 

/ /(A(xo, r))rfr 

s-so J so 



= Ai+A2-A 



3- 



Term by term: For Ai first use s — Sq = (^(s) —^{so))t/{x — y). By the assumption 
s G T and Lemma 3.4(i), we can replace s by an average over [so,s] fl T at the 
expense of an 0(s — sq) error term. Since T covers almost all of [sq, s] we can ignore 
the restriction to T in the average. These steps give 



Ai 





1 




My, 


s)dy 
1 


s 

X 


-So 

-y 



dr I X{y,T)dy + 0{s- So) 



t {s-So){^{s)-^{so)) Jso J^iso) 

— ' —^HC{so),so) 
as s \ So along points of T, by the assumption (^(sq), sq) G Leb(A). Next, 

A2^/(A(xo,so))-/(A(e(so),so)) 

as s \ So along points of T, by the assumptions sq G T and (xq, sq), (^(so), so) G 
Leb(/ o A) and by Lemma 3.4(ii). Finally, 

-43 — ^ /(A(xo,so)) 

by the assumption sq G Leb{/(A(xo, •))}• Altogether we have proved (3.16) for a.e. 

So- 

By (2.8), (3.16) implies 7'(s) > —g{{x — y)/t) for a.e. s. Since v{y,Q) = uo{y), 
we can rewrite (3.15) as 

v(x,t) ^uoiy) + / l'{s)ds>uo{y)-tQ{{x-y)/t). 



22 

This argument is valid for a.e. y, and consequently by the Hopf-Lax formula (2.10), 
v{x,t) > u{x,t). We have proved (3.14), and thereby (3.6) for x = xq. The choice 
of Xq was such that it covers a.e. point in R. 
For the converse, suppose 

(3.17) / f{X{x,s))ds= / f{p{x,s))ds for a.e. (x, t). 

Jo Jo 

(3.8) is valid for both A and p for a.e. t. Consider such a t for which also (3.17) 
holds for a.e. x. The right-hand side of (3.8) is J-^dx'i/j{x) j^ f{X{x, s))ds, so by 
(3.17) its value is not changed by replacing A with p. We conclude that for a.e. t, 
f^X{Xjt)i/j{x)dx = f^p{x,t)'i/j{x)dx for all test functions i/j. This implies X = p 
a.e. 



4. Second class particles in the variational coupling 

4.1. Construction and variational coupling. First we construct the process. 
Let {Di : i G Z} be a collection of mutually independent rate 1 Poisson jump time 
processes on the time line (0, cxd), defined on some probability space (O, JF, P). We 
construct the process r]{t) — {r]{i, t) : i G Z) on O in terms of "current particles." 
These form a process z{t) = {z{i,t) : z G Z) of labeled particles that move on Z. 
The location of the ith particle at time t is z{i,t), and these satisfy 

(4.1) 0<z{i + l,t)-z{i,t)<K. 

In the graphical construction, z{i) attempts to jump 1 step to the left at epochs of 
Di. If the jump violates (4.1), it is suppressed. We can summarize the jump rule 
like this: 

If t is an epoch of Di, then 
(4.2) 

z{i, t) = max{^(i, t—) — 1, z{i — 1, t— ), z{i + 1, t—) — K}. 

To construct z{-) from rule (4.2) we exclude an expectional null set of realiza- 
tions of {Di}. For the remainder of this section, fix a realization {Di} with these 
properties: 

(4.3a) There are no simultaneous jump attempts. 

(4.3b) Each Di has only finitely may epochs in every bounded time interval. 
(4.3c) Given any ti > 0, there are arbitrarily faraway indices io « << zi such 
that Dig and D^-^ have no epochs in the time interval [0,ti]. 

These properties are satisfied almost surely, so the evolution z{t), < t < oo, is 
well-defined almost everywhere on O. Then the process r]{t) is defined by 

(4.4) r]{i,t) = z{i,t)-z{i-l,t). 

It should be clear that r]{-) satisfies the rules described in Section 2. 

Define a family {w^ : /c G Z} of auxiliary processes on O. Each w^{t) = {w^{i,t) : 
i G Z) is a process of the same type as z{t), so (4.1) is in force for all w^{i,t). The 
initial configuration of w'^ depends on the initial position z{k, 0): 



w\i,Q) = [ 



z{k,0), i>0 



/ 7 i^\ 
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The process w'^ reads its jump commands from the same Poisson processes {-Di}, 
but with a translated index: at epochs t of Di-^-k, 

w^{i,t) = m8ix{w''{i,t-) - l,w''{i - l,t-),w''{i + l,t-) - K}. 

We prove a number of pathwise properties for the processes z and {w'^} that 
are vahd under the fixed reahzation {D^} that satisfies (4.3a-c). Since (4.3a-c) 
are vahd almost surely, the properties we derive are almost sure properties, but we 
leave the modifier "almost surely" out of the statements that follow. 

The coupling of the processes z and {ty^} preserves ordering. One consequence 
is the following inequality, valid for all /c < /, all i and all t: 

(4.5) w^{i -l,t)< w''{i -k,t) + [z{l, 0) - z{k, 0)]. 

The usefulness of the family of processes {w'^} lies in this "variational coupling" 
lemma: 

Lemma 4.1. For all i E Z and t > 0, 

(4.6) z{i,t) = snpw''{i-k,t). 

kez 

In the case K — oo the supremum in (4.6) is over {k : k < i} because w^{i—k, t) = 
— oo for k > i. 

For some purposes it is convenient to modify the processes w^ so that they start 
at level zero and advance in the positive direction. So define processes {^^} by 

eiht)=z{k,o)-w'{i,t). 

The process $,^ does not depend on 2(/c,0), and depends on the superscript k only 
through a translation of the i-index of the Poisson processes {Di}. Initially 

t. ^ r 0, i>0 

[ —Ki , z < . 

Dynamically, at epochs t of Di-^-k, 

eit^ t) = min{^'=(i, t-) + 1, ei^ - 1, t-), e{^ + 1, t-) + K}. 

These jumps preserve the inequalities 

C'"{i,t)<C''{i-l,t) and C''{i,t)<CHi + '^,t) + K. 
Now (4.6) can be written as 

(4.7) z{i, t) = snp{z{k, 0) - ^''(i - k, t)}. 

kez 

This variational formula is a microscopic analogue of the Hopf-Lax formula (2.10), 
and the basis for the proof of the hydrodynamic limit. To prove Theorem 1, one 
shows that the right-hand side of (4.7) converges to the right-hand side of (2.10). 
For more details about the construction, the variational coupling lemma, and the 

e„ ^-^ „J — .„„,. „ — „,,ij- c ■^^■^: r-innn\ 
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4.2. The second class particle. We investigate the position X{t) of a second 
class particle in the variational coupling. Fix again a sample point so that we can 
regard {-Di}, (?7(i,0)), and X{0) as deterministic, with (4.3a-c) satisfied by {Di}. 

As defined in Section 2.2, X(t) is the location of the unique discrepancy between 
two processes rj and fj that satisfy initially fj{X{0), 0) = rj{X{0), 0) + 1 and ?)(«, 0) = 
r]{i,0) for i ^ -^(0). In particular, the second class particle is not counted in the 
?7- variables and by the X-exclusion rule always r\{X{t\ {) <K — \. 

We define r\ and f] by (4.4), in terms of processes z and z that initially satisfy 

z{i, 0) = z{i, 0) for i < X(0) - 1 and z{i, 0) = z(z, 0) + 1 for z > X(0). 

Processes z and 5 follow the same Poisson processes {-Di}, so this is the basic 
coupling. The location X{f) is always uniquely defined as the location of the dis- 
crepancy, as shown in this lemma: 

Lemma 4.2. For all t > there is a finite location X{t) such that 

(4.8) z{i,t) = z{i,t) fori < X{t) -1 and z{i,t) = z{i,t) + I fori >X{t). 

Proof. Given any finite time horizon to, find indices ■■•<z_i<io<^i<"'' such 
that Di^ has no jump epochs in the time interval (0, to] for all m, im -^ ±oo as 
m -^ ±oo, and zq < X{0) < ii. Particles z{im) and z{im) do not move up to time 
to. In each finite segment (im-i, im) there are only finitely many jump epochs up 
to time to. Arrange these finitely many jump epochs in temporal order, and then 
prove by induction that the lemma holds after each jump time. D 

(4.8) defines X{t) in our framework. Next we give X{t) a representation that is 
a microscopic analogue of the characteristics. 

Proposition 4.1. We have the formulas 

(4.9) X{t) = inf{i e Z : z{i, t) = z{k, 0) - ^^i - k, t) for some k > X{0)} 
and 

(4.10) X{t) = 1 + sup{i e Z : 5(z, t) = z{k, 0) -^'"{i- k, t) for some k < X{0)}. 

Proof. Consider these four statements: 

(4.11) lii<X{t), then z{i,t) > z{k,0) - C'"{i - k,t) for all A; > X(0). 

(4.12) If z > X(t), then z{i,t) = z{k,0) -i^{i-k,t) for some k > X(0). 

(4.13) If i < X(t), then 5(z, t) = z{k, 0) - i^{i - /c, t) for some k < X(0). 

/A -\ A\ Tf „■ \ -v/-i-\ -i-1 ;:;/'„■ J.^ ^ ::^/i„ ^^ tk / ■ u J.^ e — „ii ;„ ^ x^^'^^ 
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These statements imply that the infimum and supremum in the formulas (4.9-10) 
are well-defined and finite, and also prove the formulas themselves. 
To contradict (4.11), suppose i < X{t) and 

(4.15) z{i,t) = z{k,0)-C''{i-k,t) 

for some k>X{0). Then by (4.8) 

z{i,t) =2{k,0)-^^{i-k,t)-l 

which contradicts the variational formula (4.7) for process z. This contradiction 
proves (4.11), and also implies that there must exist some k < X{0) such that 
(4.15) holds. By (4.8), (4.15) then turns into 

z{i,t) =z{k,Q)-i''{i-k,t), 

and this proves (4.13). 

(4.12) and (4.14) are proved by a similar argument that starts by contradicting 
(4.14). D 

The second class particle may jump either left or right, but a certain monotonicity 
can be found from the representation (4.9): 

Proposition 4.2. Let 

ki{t) = supjfc > X{0) : z{X{t),t) = w^{X{t) - k,t)} 

denote the maximal k that satisfies the requirement in formula (4.9), or infinity. 
Then ki{t) is nondecreasing as a function of time. If there are initially arbitrarily 
large indices j such that r]{j, 0) < K — 1, then ki{t) is finite for all t. 

We shall not prove Proposition 4.2 for we make no use of it in this paper. 

5. Proofs of the laws of large numbers 

Now assume we are in the setting described in the paragraph preceding Theorem 
1. The initial particle configurations (?7^(i,0) : i G Z), n G N, and the Poisson 
jump time processes {-Di} are defined on a probability space (0,JF, P). Define 
initial configurations (ztj,(z,0) : z G Z) by setting 2tj,(0,0) = and using (4.4). 
Construct the processes Zn{-) as in Section 4.1, and then define the processes ?7n(-) 
by (4.4) in terms of Zn{-). The location of the second class particle is given by the 
variational representation in Proposition 4.1: for the nth process 

Xn{t) = inf{z G Z : Zn{i,t) = w^{i - k,t) for some k > ^^(0)} 
= inf{z G Z : ^^(z, t) = Zn{k, 0) -^''{i- k, t) 
for some k > Xn{0)}. 

Theorems 3 and 4 rely on the basic hydrodynamic limits, which we summarize 
here from Seppalainen (1999) in terms of the ^-variables. The following statements 
all hold almost surely: 

('5.1') lim —z„(\nii].()) — ii,n(ii) for all ii G R.: 
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(5.2) lim -^['"'^'([nx], nt) = tg{x/t) for all a;, r G R and t > 0; 
and 

(5.3) lim —Zn{\nx\,nt) = u{x,t) for all x G R and t > 0. 

The function uq is defined by (2.7) from the given initial profile po- (5.1) is a 
restatement of assumption (2.3). (5.2) is proved by subadditivity, and this limit 
defines g. Then u{x,t) is defined in terms of uq and g by the Hopf-Lax formula 
(2.10). Next (5.3) is proved, by (5.1), (5.2), and the variational coupling (4.7). 
Lastly, / is defined as the conjugate of g, and ti(x, t) is identified as the viscosity 
solution of (2.11). 

The following estimate [Seppalainen (1999), Lemma 6.1] reduces the range of 
indices needed in the variational formula. It will be used several times in the 
proofs. 

Lemma 5.1. Define 

Cn{iJ,t)= max {zn{k,0)-^''{i-k,t)}. 
i-i<k<i+r 

Fix r > t. Then there exists a finite positive constant C = C{r, t) such that for all 
i &7j and all n>l, 

(5.4) P{zr,{i, nt) ^ Cu{h nr, nt)) < e"^'^. 
First we prove Theorem 3, so recall assumption (2.20): 

(5.5) lim -Xn(0) = h a.s. 

n^oo 77, 

Proposition 5.1. Assume (2.3) and (2.20). Then 

limsupn~"'^X^(?it) < x'^{h,t) almost surely. 

n — >oo 

Proof. Fix xi > x~^{bjt). We shall show that lim sup^^ ^n~^Xn{nt) < xi almost 
surely. By definition x+(6, t) = sup{a; : y~{x,t) < 6}, so y~{xi,t) > b. Choose ci 
so that 6 < ci < y~{xi,t). 
Let r be the event on which 

(5.6) Xn{nt) > nxi for infinitely many n, 

assumptions (4.3a-c) of Section 4.1 hold so that we can use variational formulas 

(4.7) and (4.9), the laws of large numbers (5.1)-(5.3) and (5.5) hold, and 

(5.7) Zn{[nxi], nt) = w^{[nxi] — /c, nt) 

for some k > [nco] for all large enough n, for some fixed finite cq. Except for (5.6) 
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Suppose r is not empty. Fix a sample point in F, and suppose n is such that 
Xn{nt) > nxi. Then Zn{[nxi],nt) ^ w'^dnxi] — k,nt) for all k > Xn{0), so there 
must exist a, k < Xn{0) such that (5.7) holds. Since limit (5.5) holds on F, X„(0) < 
[nci] for large enough n. Thus for a subsequence of n's there exists a /c in (5.7) 
such that [nco] <k< [nci]. 

Pick a partition cq = yo < ?/i < • • • < Vm = ci oi [co,ci] with mesh Ay = 
max(yi_|_i ~yi). Fix i so that for infinitely many n in the subsequence picked above, 
some k in (5.7) satisfies [nyi] < k < [nyi-^-i]. By (4.5) we have, for infinitely many 
n, 

2:^1 ([nxi], nt) = w^{[nxi] — k, nt) 
(5.8) 

< «;i7^^H[na:i] - [ny,],nt) + {zr,{[ny^+i],0) - Zn{[ny^],0)}. 

Divide by n and let n ^ oo along the appropriate subsequence of n's. Since the 
laws of large numbers hold for our fixed sample point, we get 

7 — -j + {^oiVi+i) -uo{y^)}. 

Take the mesh Ay -^ 0. By continuity and compactness, we conclude that the Hopf- 
Lax formula (2.10) for u{xi,t) has a maximizer in [co,ci], so that y~{xi^t) < ci. 
This contradiction with the choice of ci implies that F is empty. D 

Proposition 5.2. Assume (2.3) and (2.20). Then 

liminf n~^X^(?it) > x~{b,t) almost surely. 



Proof. Fix xi < x~{b, t). Consider the event F on which Xn{nt) < nxi for infinitely 
many n, assumptions (4.3a-c) hold, and the laws of large numbers (5.1)-(5.3), and 
(5.5) hold. We assume that F has positive probability, and derive a contradiction. 
(4.12) and Xn{nt) < nxi imply that for some k > X^(0), 

(5.10) Zni[nxi],nt) = Zn{k, 0) - ^^{[nxi] - /c, nt). 

As in Section 4.2, we have the processes 5^, coupled with Zn through the Poisson 
processes {-Di}, that satisfy 

Zn{i, nt) = Zn{i, nt) for i < Xn{nt) 
and Zn{i, nt) = z^ii, nt) + 1 for z > Xn{nt). 

(4.10) implies that 

(5.12) z^{\nxilnt) = z^{k, 0) - i^{\nxi] - /c, nt). 

By (4.14), (5.12) cannot hold for any k < X„(0). By applying Lemma 5.1 to the 

processes Zn and by shrinking F by no more than a null set, we may assume that 

for large enough n the k in (5.12) must lie in the range X^(0) < k < [nci] for a 
rpi u,. re: c;^ . .i,,j„ 4-u„4- e — i. — , r — . i ^ l ^ r„„ i „.i „ 
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can be taken arbitrarily close to b. By left-continuity of x~{-, t) and the assumption 
xi < x~{h,t), we may choose cq so that xi < x~{co,t). 

Now starting from (5.12) repeat the argument around (5.8)-(5.9) for Zn- By 
(5.11) Zn satisfies the limits (5.1) and (5.3). As in the step after (5.9), we get that 
u{xijt) has a maximizer in [co,ci] so that y^{xi,t) > cq. By definition (2.18) this 
implies x~(co, t) < xi and contradicts the choice of cq. □ 

Propositions 5.1-5.2 prove Theorem 3. 

Proof of Theorem 4(Q')- Assume x{y, t) = x^{y, t) for all y eH. Let F be the event 
on which (4.3a-c) and (5.1)~(5.3) hold, and also for large enough n 

(5.13) -nbo < Xr,{0) < nbo 

for some constant Bq. We prove that 



(5.14) lim 



>CXD 



Xnjnt) ^ / X^(0) ^ 



0. 



n \ n 

almost surely on F. By assumption (2.23), Theorem 4(a) follows by repeating this 
for a sequence of 6o's increasing to oo. 

Since x~{y,t) is left-continuous and x^{y,t) right-continuous in y, the assumed 
unique characteristic x{y,t) = x'^{y,t) is continuous in y. Given e > 0, choose a 
partition — 6o = ro < ri < ■ ■ ■ < r/^ = 6o so that x{rjj^i,t) < x{rj,t) + e/2. Set 
Xj — x{rjjt) +e. Shrink F by a null set to guarantee that a.s. for every j = 0, . . . , /i, 
for large enough n there exists a, k > nco such that 

Zn{[nxj],nt) = w^{[nxj] — k^nt). 

This can be done by Lemma 5.1 if cq is chosen small enough. 

To contradict (5.14), suppose F contains a sample point for which 

(5.15) Xn{nt) > nx{n~^Xn{0), t) + ne 

for infinitely many n. Restrict this subsequence of n's further so that for some fixed 
h 

(5.16) nrj < X„(0) < nr^+i 

for all these n. By the monotonicity of a:(-,t), (5.15-16) imply that Xn{ni) > nxj 
for infinitely many n. For each such n there exists k < Xn{0) and no k > Xn{0) 
such that Zn{[nxj], nt) = w^([nxj] — k, nt). By the assumptions made above, ucq < 
k < nrj+i. Repeat the partitioning and limiting argument around (5.8)-(5.9) to 
conclude that y~{xj,t) < rj+i, and then by definition x"'"(rj+i, t) > Xj. Now we 
have contradicted Xj = x{rj, t) + e > a;(rj+i, t) + e/2. 

It remains to show that assuming X„(nt) < na:(n~^X^(0), t) — ne for infinitely 
many n also leads to a contradiction. Combine this with (5.16) to get 

Xn{ni) < nx{n~^Xn{0),t) — ne < na:(rj_|_i, t) — ne = nx'j. 

As in the proof of Proposition 5.2, switch to the Zn processes to get 

Zn{[nx'j],nt) = Zn{k,0) -^''([nx^] -k,nt) 

for some k such that nrj < ^^(0) < k < nci, where ci is a fixed constant. Prom this 
conclude that y~^{x'j, t) > rj, which implies x{rj, t) = x~{rj^ t) < x'^ = x(rj_(_i, t)— £, 
a contradiction. D 
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Lemma 5.2. Let c be a finite constant. 

(a) Suppose that Xn{0) > nc for large enough n, almost surely. Then for any 
X < c — t, Xn{nt) > nx for large enough n, almost surely. 

(b) Suppose that X„(0) < nc for large enough n, almost surely. Then for any 
X > c + t, Xn{nt) < nx for large enough n, almost surely. 

Proof, (a) Pick £ > so that x < c — t — 2e, and set ci = c — e. 

As an intermediate claim, we prove that almost surely, for large enough n, 

(5.17) ^''(M - /c, nt) = -K ■ {[nx] - k) for all k > [nx] + [n{t + e)]. 

In other words, the claim is that ^^{[nx] — /c, ■) has not jumped by time nt. The 
time of the first jump of ^'^{[nx] — k,-) is distributed as X]i=o ^' where {Yj} 
are i.i.d. rate 1 exponential random variables. The estimate Pi^^^iYj < Ns) < 
exp[— A^k(s)] is valid for < s < 1. The rate function is k{s) — s — 1 — logs for 
s > 0, which satisfies k' {s) < < k{s) for < s < 1. Thus 

(k — [nx\ 
Y^ Yj < nt 
j=o 

< J2 exp[-m(nt/i)] <Cie-^2'" 

i>[n{t+e)] 

for finite positive constants Ci, C2. This proves the claim, and (5.17) holds for large 
enough n a.s. 

Since r)n{X{0), 0) < K - 1 and [nci] < X(0), 

k 

Zn{k,0)-Zn{[nci],0)= Yl Vn{i,0)<K{k-[nci])-l 

i=[nci] + l 

for any k > X„(0). For the fixed x any k > Xn{0) falls within the range in (5.17) 
for large enough n. So for large enough n, 

Zn{k, 0) - ^''{[nx] - k, nt) = Zn{k, 0) + K{[nx] - k) 
< z^{[nci],0) + K{[nx] - [nci]) - 1 = Zn{[nci],0) - ^["^^^(M - [nci],nt) - 1. 

The last equality comes from (5.17) again because [nci] also falls within the range of 
/c's in (5.17) for large enough n. The string of inequalities shows that no k > ^^(0) 
can give the supremum in the microscopic variational formula (4.7) for ^^([^^2;], nt), 
because each k > X„(0) is inferior to [nci]. But if Xn{nt) < [nx], by (4.12) some 
k > Xn{0) would be a maximizer. This contradiction shows that Xn{nt) > [nx] for 
large enough n, a.s. 

The proof of part (b) is similar, but instead of (5.17) use that ^'^([nx] — k, nt) = 
for all k < [nx] - [n{t + e)]. D 

Proof of Theorem 4(b)- The goal is to show that 



(5.18) lim !^^^-J^^^;s,t) 



= 
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Let uq{x) = u{x,s), u{x,t) = u{x, s + t) for r > 0, and let 3;(6; 0, r) denote 
characteristics for the solution u. Then x^{b;s,t) = 3:^(6; 0,t — s). From the 
assumptions and Proposition 3.1 we know u has unique characteristics: x{b; 0, r) = 
x±(6;0,r). 

Let Zn{ijt) = Zn{i,ns + t) be the nth process restarted at time ns. Let Xn{t) 
be the location of a second class particle in the process z^, with initial location 
-^n(O) = Xn{ns). In restarted terms, (5.18) is the same as 



(5.19) lim 



Xn(n(t-.))__(M!l^,0,t-s 



n \ n 



a.s. 



(5.19) follows from Theorem 4(a) if we check that the ingredients of the proof are 
valid for the processes Zn and Xn- 

(i) We can obtain the variational formulas (4.7) and (4.9) for Zn and Xn in a form 
where the ^-part depends on n through the time shift. For any fixed restarting time 
T > 0, let 9t- denote the restarting operation on the probability space of the Poisson 
processes {Di}: if the epochs of Di in (0, oo) are ti < t2 < ts < ■ ■ ■ with tj-i < 
T < tj, then the epochs of 9rDi in (0, oo) are tj — t < tj+i — t < tj+2 — t < ■ ■ ■ . 

The same proofs that originally gave Lemma 4.1 and Proposition 4.1 also yield, 
for t > T, 

(5.20) z{i, t) = sup{z{k, r) - C'^ii - k,t - r) o 9^} 

k 

and 

(5.21) X{t) = inf{i : z{i, t) = z{k, r) - ^^{i - k,t - r) o 6^ for some k > X{t)}. 
Now for process Zn take r = ns, and (5.20)-(5.21) give for t > 

Zn{i, t) = SVip{Zn{k, 0) - ^'^ {i - /C, t) O Ons} 
k 

and 

Xn{t) = inf{z : Zn{i, t) = Zn{k, 0) — ^ {i — k,t) o 0ns for some k > Xn{ns)}. 



(ii) Limits (5.1) and (5.3) are true by definition of Zn and u. Limits (5.2) for 
^[nr]Qj^^j^ jil^oOns are true by the proofs in Seppalainen (1999) because these limits 
are derived from summable deviation estimates and the Borel-Cantelli lemma. The 
distribution of ^["^^1 is not changed by the time shift, so the same deviation estimates 
are valid. 

(iii) Finally assumption (2.25) for X^(0) = Xnijis) is a consequence of (2.25) for 
X^(0) and Lemma 5.2. D 
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